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Abstract 

We present a parametric family of Riordan arrays, which are obtained by multiply¬ 
ing any Riordan array with a generalized Pascal array. In particular, we focus on some 
interesting properties of one-parameter Catalan triangles. We obtain several combina¬ 
torial identities that involve two special Catalan matrices, the Chebyshev polynomials 
of the second kind, some periodic sequences, and the Fibonacci numbers. 


1 Introduction 

There are two infinite lower triangular matrices in the mathematical literature, that are 
both called Catalan triangles. Let us denote them by C and B. respectively. By way of 

1 Corresponding author. 

2 Work performed within the activities of Centro de Analise Funcional, Estruturas Lineares e Aplicagoes 
(Faculdade de Ciencias, Universidade de Lisboa) and supported by the fellowship SFRH/BPD/48223/2008 
provided by the Portuguese Foundation for Science and Technology (FCT). 


1 



illustration, their first rows are shown below: 
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Note that the 0 th column of C is given by the sequence of Catalan numbers (C n ) n > o, 
where C n = yy-j-( 2 ”). In turn, the 0 th column of B is given by (C n+ i) n >o- This similarity 
explains why they are called Catalan triangles. 

According to the online encyclopedia of integer sequences (OEIS) [11], the non-zero 
entries of C (sequence A033184) are given by the formula C n>k = |Ty ( 2n ~ k ) , for n > k > 0. 
They are known as the (ordinary) ballot numbers. Aigner [2] used these numbers (and 
generalizations of them) to enumerate various combinatorial instances. Furthermore, Ferrari 
and Pinzani [8] gave an interpretation of C , using the ECO method, and a suitable change 
of basis in the vector space of one-variable polynomials. On the other hand, the non-zero 
entries of B (A039598) are given by the formula B nk = “yy Cn-k) j for n > k > 0. The 
numbers B nk first appeared in the work of Shapiro [17, 19], in problems connected with 
non-intersecting paths, random walks, Eulerian numbers, runs, slides, and moments. 3 We 
will refer from now on to C as the Aigner array, and to B as the Shapiro array. 

The aim of this paper is to discuss some interesting features regarding the product matrix 
R(r) = RP(r). where R is a given Riordan array, and the parameter r is any real or complex 
number. The notation P{r) stands for a generalized Pascal array, such that -P(O) = I is the 
identity matrix, and P (1) = P is the classical Pascal array of binomial numbers, that comes 
arranged as a lower triangular matrix. We elaborate, in particular, on the product matrix 
C(r ) = CP{r). Clearly, we have C(0) = C. It turns out that C(l) = B. The factorization 
B = CP resembles Barry’s definition [3, see Sections 6.3 and 6.5] of the generalized ballot 
array given by Bal = Cat • Bin; where Bin = P , and Cat = diag(l, C ). 

Recently, Yang [23] introduced a generalized Catalan matrix, given by 


C[a, 6; r] = 


1 — y/1 — 4 rz\ a (1 — \/l — 4rz 


2 rz 


2 rz 


where a, b are integer numbers, and r is an arbitrary parameter. Such a generalization 
provides a unified way of presenting the matrices C and B, along with many other matrices 
related to Catalan numbers. One can check that the 0 th column of C[a, 6; r] forms a sequence 

3 Actually, Shapiro used the formula = ^(A), for n,k> 1. His formula is just a shifted version of 
the formula shown previously for B n ,k- 
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of degree n monomials in r, whose coefficients are polynomials in a and b, that have rational 
coefficients. The one-parameter Catalan triangles that we will study are instead given by 

/1 — 2 rz — y/l — 4z 1 — 2 rz — y/1 — 4z\ 

\ 2z(l — r + r 2 z) ’ 2(1 — r + r 2 z) ) 

We shall show explicitly that the 0 th column of C(r ) forms a sequence (C'(r) n 0 ) ri > 0 of degree 
n monic polynomials in r. These polynomials have positive integer coefficients. In particular, 
the constant term of C(r) n o is the Catalan number C n . For this reason, we shall refer to 
C(r) as a Catalan triangle too. 

It is also worth mentioning that He [9] introduced recently, a family of Catalan triangles 
which depend on two parameters. His (c, r)-Catalan triangles are based on the sequence 
characterization of Bell-type Riordan arrays. We emphasize that they are also different from 
the Catalan triangles C(r). One way to see this difference is to compare, once again, the 0 th 
column of both types of matrices. On the one hand, one can check that the 0 th column of 
a (c, r)-Catalan triangle forms a sequence of degree n homogeneous polynomials in c and r, 
whose coefficients are the well-known Narayana numbers. On the other hand, the 0 th column 
of a Catalan triangle C (r) forms a sequence of one-variable polynomials of degree n, whose 
coefficients are the ordinary ballot numbers. Recall that both the Narayana and the ballot 
numbers are known to refine the Catalan numbers. 

The fundamental theorem of Riordan arrays, and the isomorphism between Riordan 
arrays and Sheffer polynomials (the former being the coefficient matrices for the latter), 
allow us to interpret C'(r) as a change of basis matrix in the space of polynomials in one 
variable. Thus, for instance, the Shapiro array B is seen as the change of basis between the 
Chebyshev polynomials of the second kind, and the standard polynomial sequence of powers 
of one variable. In this way, by specifying values for the one-variable of the corresponding 
polynomials, several formulas related to periodic sequences, the sequence of natural numbers, 
and the sequence of Fibonacci numbers, can be stated and generalized. 

2 General setting 

The literature on Riordan arrays is vast and constantly growing (see, for instance, [1, 4, 
5, 7, 9, 10, 12, 13, 14, 15, 16, 18, 20, 21, 22]). In this section we recall the fundamental 
theorem of Riordan arrays, and discuss some general facts regarding the generalized Sheffer 
polynomial sequence associated with any Riordan array. Then we define a one-parameter 
family of Riordan arrays, and focus our attention on Catalan arrays in Section 3. 

2.1 Riordan Arrays 

Let d(z ) = d 0 + diz + d 2 z 2 + • • •, and let h(z) = h 0 + h ± z + h 2 z 2 + • • • be two formal 
power series in z with coefficients in a given integral domain R , with unit 1. Assume that 
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do = d(0) 7 ^ 0, and h 0 = h( 0) = 0. A Riordan array R is an infinite lower triangular matrix, 
whose entries are given by 

Rn,k = [Z n }d{z)h{z) k . (1) 

The operator [z n ] acts on a formal power series f(z) = J2 k>0 fkZ k by extracting its n th 
coefficient; that is, we have [ z n ]f(z ) = f n . It is customary to denote R = ( d(z),h(z )). 
In addition, the set of Riordan arrays for which h'{ 0) 7 ^ 0, and that is equipped with the 
multiplication given by 

(di(z), h 1 (z))(d 2 (z), h 2 {z )) = (d 1 (z)d 2 (h 1 (z)), h 2 (h 1 (z))), ( 2 ) 


forms a group. We denote this group by fHio. The group identity is the usual identity matrix 
I = (1 ,z). Furthermore, for any R G 94io, its inverse is given by R = (1 /d(h(z)),h(z)), 
where h(z) denotes the compositional inverse of h(z). 

The order of a formal power series g(z) = Y2 k> o 9kZ k is the minimal index k such that 
g k 7 ^ 0. Since i?[[z]] denotes the ring of formal power series in z with coefficients in R , 
we will write R r [[^]] to denote the set of formal power series of order r. Therefore, if 
R = ( d(z),h(z )) G iHio, then d(z) G i?o[[z]], an d h(z) G [[*]]. For simplicity, and without 
loss of generality, we further assume that do = 1 , and hi = h'{ 0 ) = 1 , so that R n . n — hi — 1 , 
for all n > 0 . 

Now, let us define an action * of the group fHio on i?o[[^]] by 

(d(z), h{z)) * g(z) = d(z)g(h(z)). (3) 


Formula (3) represents the traditional matrix-column multiplication. Therefore, it is also 
equivalent to saying that, for any n > 0, the following identity holds 

n 

R n,k9k = [z n }d(z)g(h(z)). (4) 

k=0 

Equation (4) is known as the fundamental theorem of Riordan arrays. It is usually used to 
deal with sums and combinatorial identities. In particular, let g(z) = 1 + xz + x 2 z 2 + • • ■ = 
. We write 

1 —XZ 

= = . (5) 

k =0 ^ ' 

We call (Pn(x)) n > o the Sheffer polynomial sequence associated with R. Formula (5) implies 
that the generating function of Pn{ x ) is given by 

, i(Z l v (6) 

s 1 - xh ^ 


Likewise, the generating function for the Sheffer polynomial associated with the inverse array 
R = (l/d(h(z)), h(z)) is 


EpF'm*” 

n> 0 


1 1 

d{h{z)) 1 — xh(z) 


(7) 
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Since R is invertible, and the sequence (x n ) n > 0 is a basis for the linear space C[x], the 
sequence {Pni x ))n >o is also a basis for C[x]. Therefore, any R e 9tio can be interpreted as 
the change of basis matrix between the standard basis (x n ) n >o, and the basis of polynomials 
(Pn( x ))n> o- Hence, we have 

n 

Y R ^Pk o) = x ™■ ( 8 ) 

k=0 

We will use Equation ( 8 ) in Section 4 to deduce some interesting combinatorial identities 
involving the Riordan arrays C and B. 


2.2 One-parameter Riordan arrays 

The classical Pascal array P of binomial numbers (”) is a Riordan array, since we can 
write P = (p(z), zp(z)), where p(z) = U—. Given any real number r, the Riordan array 
P(r) = (p(rz), zp(rz)) = (yU—, jfqx) is known as the generalized Pascal array of parameter 
r. Its entries are given by P(r) n ^ = \£)r n ~ k . 

Definition 1. Let R be in 9tio. The r-Riordan array R{r) is defined by R(r) = RP{r). 

Set R = (, d{z),h(z )). Then, by Equation (2), we have R(r) = , rrlpy) • Using 

Equation (1), it is easy to check that the entries of R(r) are given by 


R(r)n,k = [-1 


d(z) 


h(z) 


1 — rh(z) \ 1 — rh(z) 


= Rt 


d(z)h(z) h 


(l — rh(z)) k+1 

Lemma 2. For any real number r, the Sheffer polynomial sequence (p^ r \x )) n >o satisfies 
the identity 

Pn {r \ x ) =Pn( x + r)- 

Proof. By Equation (5), we have 

d{z) 

«r) n , t x* = [V] 1 ' = [ z «] _ 


Pn {r \ X ) =Y R ( 


= p R (x + r). 


k=0 


1 — X 


h (z) 

1 — rh(z) 


(x + r)h{z ) 


□ 


Corollary 3. Let n be a non-negative integer. The constant term and the sum of coefficients 
of Pn ir) {x) are given by 


Pn ,r \ 0 ) = R(r )„,o -Vi? 




(9) 


k=0 


p\ 


R ^(l)=Y R ( r )n,k = R(r+l)n, 0 . 


( 10 ) 


k=0 
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3 r-Catalan triangles 


Let us recall that the generating function for the well-known Catalan numbers C n = ( 2 (‘) 

is c(z) = We apply Definition 1 to the Riordan array C = (c(z), zc{zf). 

Definition 4. The r-Catalan triangle C(r ) is given by C(r) = CP(r). Then 


C(r) 



c(z) 

— rzc(z )’ 


zc{z) \ 
1 — rzc(z)J 


f 1 — 2 rz — y/l — Az 1 — 2 rz — y/l — 4 z\ 
\ 2z(l — r + r 2 z) ’ 2(1 — r + r 2 z) ) 


( 11 ) 


In particular, we have (7(0) = (7. Moreover, since the identity c(z) = 1 + zc(z) 2 holds, 
it follows that (7(1) = CP(1) = (j^j, = (c(z) 2 , zc(z) 2 ) = B. 

The next formula follows from Definition 4. 


Lemma 5. The entries of C{r ) are given by 


C' {r)n,k 



c(z) 

1 — rzc(z) 


fe+i 




i + k + 1 

n+1 


^ 2 n — i — 



( 12 ) 


The matrix C(r) is an array of polynomials in r. By way of illustration, the first entries 
of C(r) are 


C(r) 


( 1 

1 + r 1 

2 + 2r + r 2 2 + 2r 1 

5 + 5r + 3r 2 + r 3 5 + 6r + 3r 2 3 + 3r 1 

14 + 14r + 9r 2 + 4r 3 + r 4 14 + 18r + 12r 2 + 4r 3 9 + 12r + 6r 2 4 + 4r 1 


\ 


V 




For every k > 1, column k of C{r ) is the (k + l)-fold convolution of column 0. Thus, for 
example, the first columns of (7(2) are given by the following sequences: 


(7(2) = ( A001700 A008549 A045720 A045894 R035330 • • •) . 


The 0 th column of C{r) forms a sequence ((7(r) n 0 ) n >o of monic polynomials, which are 
given by the matrix-column product (7[l,r, r 2 ,.. ,] T , as observed in Equation (9). We can 
obtain some known sequences by assigning values to the parameter r. For instance, we have 
(C(3) ri , 0 )„>o = 21049027 and (C(4) n0 )n>o = 21076025 . Observe that, by Equation (12), 
the constant term of each polynomial in the sequence ((7(r) n! o) n >o is equal to the Catalan 
number C n . Moreover, we also have 


C n+ i — 


2n + 2 


n + 2 \ n + 1 


n + 1 


n 


= B 


n,0 


= (7(l) n ,0 = 


“ n+ 1 

2=0 


n 


k =0 
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This is a well-known result. In this sense, we say that the polynomial sequence (C(r) n 0 ) n > 0 
is a refinement of the shifted Catalan sequence (C n+ i) n >o. 

Now, set h(z) = zc(z). We have z = h(z)c(h(z)). It is easy to check that h(z) = z — z 2 . 
Hence, we have C = (Yc(zc(z)), zc(z)) = (1 — z, z — z 2 ). Since C(r ) = P(r ) C , it 
follows that 


C(r)~ = 


;) I 1 - *-■ z - + 


1 + rz' 1 + rz / 

Lemma 6. The entries of C(r) 1 are given by 


f 1 + (r — 1 )z z + (r — l)z 2 \ 
\ (1 + rz) 2 ’ (1 + r^) 2 ) 


(13) 


C(r) 


-1 

n,k 



fl + (r — l)z 
\ (1 + rz) 2 


fc+i 




(14) 


Lemma 5 and Lemma 6 show that both C(r) n ^ and (— l) n ~ k C{r)^ k are polynomials in 
r with positive integer coefficients. In addition, note that the polynomial C(r) n+11 is of 
degree n, has constant term C n+ 1 , and its leading coefficient is n + 1. Both lemmas are 
useful to deduce several formulas and combinatorial identities. For instance, letting r = 0 in 
Equation (14), we obtain 




Similarly, letting r = 1 in Equation (14), we get 



Thus, the following formula is immediate. 

Lemma 7. 



The next non-trivial combinatorial identity follows directly from Equation (10) and Equa¬ 
tion ( 12 ). 

Lemma 8. 


n n—k 

EE 

k =0 i =0 


i + k + 1 (2n 

n+1 


n 


i + k 
k 


r = 


E 

i=0 


i + 1 
n + 1 


2 n — i 
n 


1 )\ 
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As for the generating functions for the Sheffer polynomials associated with C(r) and 
C(r) , the next formulas are straightforward consequences of Equation ( 6 ) and Equation 
(7), using Formula (11) and Formula (13), respectively. 

Proposition 9. 

J2Pn ir) ( X ) Z '‘ 

n> 0 


c{z) 

1 — (x + r)zc(z) 


X>? (rr, 0)z” = 


1 + (r — 1 )z 


n>0 


1 + (2 r — x)z + (r 2 — (r — 1 )x)z 2 


(16) 


In Section 4, we will use Proposition 9 to obtain more combinatorial identities. 


4 Special identities 

4.1 Identities involving the Aigner array 

Set r = 0 and x — — 1 in Equation (16). We obtain 

1 + I>f (-!)(-*)" = ++j 

n> 1 

= 1 + 2z + 3z 2 + 5P 3 + 8 z 4 + I3z 5 + 21z 6 + 34z 7 + 55z 8 + • • • 


= F 2 + F 3 z + F 4 z 2 + F 5 z 3 + ■■■ , 


where F n is the n th Fibonacci number, for n > 2. Hence, for n > 0, it holds that 

Since pf !>: ) = V"_,i we have 



(l\ 

-1 


f f 2 \ 

-F's 


( f 2 \ 

-f 3 


/ i \ 

-l 

c 1 

1 

-1 

V : ) 


f 4 

-f 5 

\ '■ / 

, or equivalently, C 

f 4 

-f 5 

V : ) 


i 

-l 

V ; / 



71+2 5 


and 


k=o 


k + 1 
n + 1 


2 n — k 
n 


-i r~ k F k+2 = i. 


Then 














4.2 Identities involving the Shapiro array 

Set r = 1 in Equation (16). Shapiro, Woan, and Getu [19] observed that the Catalan triangle 
B verified the following interesting identity: 

/ 1 \ / 1 \ 


B 


2 

3 

4 

5 


V : 7 


4 

4 2 

4 3 

4 4 

V ; 7 


(17) 


The problem of finding a matrix identity that extends the sequence (1,4, 4 2 ,4 3 ,...) to the 
sequence (1, k, k 2 , k 3 ,...) in Equation (17), was studied by Chen, Li, Shapiro and Yan [5]. 
They generalized the following known recurrence relation for B\ namely, 

^n,k -E^n—l,k— 1 Y 2 B n _\j~ T B n _\ , 

by constructing a matrix M that satisfies the identity 

( 


M 


\ 

The first column of M (that is, its 0 th column) was then interpreted in terms of weighted par¬ 
tial Motzkin paths. Here, we approach the problem of extending the sequence (1,4,4 2 ,4 3 ,...) 
to the sequence (1, k, k 2 , k 3 ,...), by keeping C( 1) = B fixed, and using Sheffer sequences. 
More explicitly, by Equation (16), we have 


1 \ 


( 1 \ 

1 + t 


k 

1 + t + t 2 

— 

k 2 

1 + t + t 2 + t 3 


k 3 



V : 7 


1 + J>n {X)z n = ~——— 


n> 1 


Furthermore, by Equation ( 8 ), we also have 


B 


/ i 

\ 


(1\ 

pf (x 



X 

pf ( x 



x 2 

\ '■ 

) 


V/ 


, or equivalently, B 


1 — (x — 2 )z + z 2 ' 


( 1 \ 

x 

x 2 

\ '■ 7 


/ i \ 

Pl_(x) 

P2 ( X ) 


V 


(18) 


7 


By Identity (15), the polynomial p^ (x) is given by 

n n 


B \ 

Pn = 


n—k 


k =0 


k =0 


n + k + 1 
n — k 


x 


(19) 


9 



















Note that pf ( x ) = where (U n (x)) n > o denotes the sequence of Chebyshev 

polynomials of the second kind [ 6 ]. Moreover, the sequence (pf (x)) n >o satisfies the three- 


term recursion 


Pn ( x ) = {x- 2 )p^_ 1 (x) -pf_ 2 (x), n> 2, 


with initial values pf (x) = 1 and pf (x) = x— 2. Whenever x G Z, the polynomial sequence 

(pf i x ))n> o reduces to an integer sequence. By specifying the values of x in Equation (18), 
we obtain many interesting relations. 

-1 -1 

Example 10 (Periodic sequences). When x — 2, the recursion pf (2) = —pf_ 2 (2), with 

initial values pf (2) = 1 and pf (2) = 0, is solved by the period-4 sequence 1,0,—1,0. 
Hence, we have 


/ 1 
2 1 
5 4 1 

14 14 6 

42 48 27 

v : : : 


\ 




( i \ 
o 

-l 

o 
i 

V ; J 


( i \ 
2 
2 2 
2 3 
2 4 

V ; 


When x = 3, the recursion pf (3) = p^_ 1 (3) — pf_ 2 ( 3), satisfying the initial conditions 
Po (3) = Pi (3) = 1, is solved by the period -6 sequence 1,1, 0, —1, —1, 0. Thus, we have 

\ 


( 1 
2 1 
5 4 1 

14 14 6 

42 48 27 

V : ; ; 


( 1 \ 
1 
o 

-i 
-i 

V ; J 


(1 \ 

3 

3 2 

3 3 

3 4 

V ; 


These two cases are better understood by recalling that U n (cosy ) = blll ^ i ^ 1 ^ . Hence, we have 

pf (2) = U n (cos -^) and pf (3) = U n (cos -|). These identities explain why we obtained 
periodic sequences. 

-1 -1 -1 

Example 11 (Natural numbers). When x = 4, the recursion p^ (4) = 2p^_ 1 (4) — p^_ 2 ( 4), 
with initial values pf (4) = 1, pf (4) = 2, is solved by the sequence (pf (4 )) n > 0 = 
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(n + l) n >o- Thus, we obtain the Identity (17): 

/ 1 \ 

2 1 
5 4 1 

14 14 6 1 
42 48 27 8 1 

V ; ; ■■■■■■} 


/1 \ 
2 

3 

4 

5 

V ; ) 


( 1 \ 
4 

4 2 

4 3 

4 4 

V ; 


The previous matrix identity is equivalent to Chen’s combinatorial formula [5]: 


£ 

k =0 


(k + l) 2 (2n + 2 


n + 1 V n — k 


— 4 n . 


Example 12 (Fibonacci numbers). When x = 5, the recursion p ^ (5) = 3p®_ 1 (5) —p^f_ 2 (5), 

with initial conditions pjf (5) = 1 and pf (5) = 3, is solved by the subsequence of Fibonacci 

numbers (p^ (5)) n >o = (i 7 2 n+ 2 )n>o- Thus, Shapiro’s triangle B , and the Fibonacci numbers 

of even index, are related by the following identity: 

/ 1 \ 

5 

5 2 

5 3 

5 4 

V ; 


Equivalently, we have 



/ 1 

\ 


/1 \ 


2 

1 


3 


5 

4 1 


8 


14 

14 6 1 


21 


42 

48 27 8 1 


55 


V ; 

; ; : : J 


V : / 

have 





k + 1 /2n + 2\ 



77 , + 1 V 77 , — k y 

^ 2 fc +2 = 


( 20 ) 


Using the well-known binomial recurrence relation (”) = + ( n fc 1 ) in Equation (19), 

and then inverting Formula (20), we obtain 


F'in+2 - (—l) n ^2 


k =0 


77 + k + 2 
77 — k 


n + + 1 

77 . — k — 1 


(-5)*. 


Moreover, seeing as F 2 n +3 = ^ 2 (n+i )+2 — E 2n+2 , we can deduce an analogous formula for the 
Fibonacci numbers of odd index; that is 


*Wi = (-i)”£ 


fc =0 


77 , + k + 2 
n — k 


77 , + k 
77 , — k — 2 


(-5) fc . 
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